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Abstract We present the Minkowski space solutions of the inhomogeneous Bethe- 
Salpeter equation for spinless particles with a ladder kernel. The off-mass shell scat- 
tering amplitude is first obtained. 
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1 Introduction 

The inhomogeneous Bethe-Salpeter (BS) equation in Minkowski space [l] 
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provides a covariant four-dimensional description of two-body scattering states. In the 
center of mass frame, and for a given incident momentum ps, the partial wave off-mass 
shell amplitude Fi depends on two scalar variables pq and |p|. It will be hereafter 
denoted by F;(po,P;Ps) setting p = |p|, Ps = |ps|. 

The on-shell amplitude F°" = Fi{pQ = 0,p — ps\Ps) determines the phase shift 
according to: 



:{i + f^Fr) (2) 



with Ep^ = ym^+pf. The knowledge of this function in the entire domain of its 
arguments - i.e. the off-shell amplitude - is mandatory for some interesting physical 
applications, like for instance computing the transition e.m. form factor 7*ci — np or 
solving the BS-Faddeev equations. This quantity has not been obtained until now. 
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The numerical solution of the BS equation in Minkowski space is complicated by 
the existence of singularities in the amplitude as well as in the integrand of These 
singularities are integrable in the mathematical sense, due to ie in the denominators 
of propagators, but their integration is a quite delicate task and requires the use of 
appropriate analytical as well as numerical methods. 

To avoid these singularities, the BS equation was first solved in Euclidean space. 
These solutions provided on-shell quantities like binding energies and phase shifts [5]. 
However we have shown 3 that the Euclidean BS amplitude cannot be used to calculate 
electromagnetic form factors, since the corresponding integral does not allow the Wick 
rotation. One therefore needs the BS amplitude in Minkowski space. 

This amplitude has been computed for a separable kernel (see [3] and references 
therein). For a field-theory kernel - the ladder and the cross ladder - it has been first 
obtained in our preceding works [5l[6] for the bound state problem. To this aim, we 
developed a method based on the Nakanishi integral representation of the BS ampli- 
tude. A similar method for the scattering states has been proposed in [Tj although the 
numerical solutions are not yet available. 

We present in this contribution a new method providing a direct solution of the 
original BS equation. It is based on a scrutinized treatment of the singularities and 
allows us to compute the corresponding off-shell scattering amplitude in Minkowski 
space. We will give the low energy parameters in the case of spinless particles and 
ladder kernel. 



2 Method 

There are four sources of singularities in the r.h.-side of the BS equation ([1} which are 
detailed below. 

(i) The constituent propagators have two poles, each of them represented as: 

— = PV—7-^ h i7r5(po — a.) 

where PV means the principal value. In the product of four pole terms, the only non 
vanishing contributions come from the product of four PVs without delta- functions, 
from the terms with three PVs and one delta-function and from the term with two 
PVs and two delta's. After partial wave decomposition the 4D integral BS equation 
is reduced to a 2D one. Integrating over pg, we obtain in addition to the 2D part, a 
ID integral over p' and a non integrated term. The singularities due to the PVs are 
eliminated by subtractions. 

(ii) The propagator of the exchanged particle has the pole singularities which, 
after partial wave decomposition, turn into logarithmic ones. Their positions are found 
analytically and the numerical integration over p'q variable is split into intervals between 
two consecutive singularities, namely: 

psingi j>sing2 psings 

[. . .] dp'o = / [...]dp'o+ [■■ ■] dpQ + [...]dp'o + --- 

Jo J singi J sing^ 

Each of these integrals is made regular with an appropriate change of variable. We 
proceed in a similar way for the p' integration. 

(Hi) The inhomogeneous (Born) term is given by the ladder kernel and is also 
singular in both variables. The positions of these singularities are analytically known. 
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(iv) The amplitude F itself has many singularities, among which the Born term 
contains the strongest ones. This makes difficult its representation on a basis of regular 
functions as well as its numerical integration. To circumvent this difficulties we made 
the replacement F = Vf, where / is a smooth function. After that, the singularities of 
the inhomogeneous term are canceled. We obtain in this way a non-singular equation 
for f which we solved bv standard methods. 




Fig. 1 On left: phase shift calculated via BS equation (solid curve) are compared to the non- 
relativistic results (dashed curve) for fi = 0.5. On right: the same comparison for the scattering 
length ao vs. the coupling constant a. 



3 Numerical results 

We first applied this method to solve the bound state problem by dropping the inho- 
mogeneous term in ([T]). The binding energies coincide, within four-digit accuracy, with 
the ones calculated in our previous work '5^ and with the Euclidean space results. 

The S-wave off-shell scattering amplitude Fq was calculated and the phase shifts 
extracted by means of eq. ([2]). They include an imaginary part, which has been also 
found, above the first inelastic threshold pi (/i) = mfi + By performing a Wick 

rotation in ((1} we derive an Euclidean space equation similar to one obtained in [5]. 
The phase shifts found by these two methods coincide with each other within 3-4 digits. 
Furthermore, the imaginary part of the phase shifts vanishes with high accuracy below 
threshold. The unitarity condition is not automatically fulfilled in our approach, but 
appears as a consequence of handling the correct solution. It thus provides a stringent 
test of the numerical method. Our results reproduce the phase shifts given in [2] within 
the accuracy allowed by extracting numerical values from published figures. 

Figure [TJ left panel, shows the phase shifts calculated via BS equation (solid curve) 
and via the Schrodinger one with the Yukawa potential (dashed curve) for the con- 
stituent mass m = I, exchange mass /i = 0.5 and coupling constant a = g'^ / (IGnm^) = 
0.5. Right panel shows the same comparison for the scattering length ao as a function of 
the coupling constant a. In the vicinity of q ~ 0.8 (for Schrodinger) and a ~ I (for BS) 
the coupling constant crosses the critical value corresponding to the appearance of a 
bound state. At this point the scattering length becomes infinite and then changes the 
sign. One can see that the differences between relativistic and non-relativistic results 
are rather large even for small and zero incident momenta. These differences increase 
with the value of a. 



Fig. 2 On left: real part of the off-shell amplitude F(po,P',Pb) for ps = 0.5, fi = 0.5. On right: 
imaginary part of F{p(),p;ps). 

We have displayed in Fig.[2]the real (left panel) and imaginary (right panel) parts of 
the ofT-shell scattering amplitude Fq{pq,p;Ps) vs. po and p calculated for ps — ^j. — 0.5. 
Its real part shows a non trivial structure with a ridge and a gap resulting from the sin- 
gularities of the inhomogeneous term. Its on-shell value F™ = -Fo(po = 0,p — ps;ps), 
determining the phase shift calculated previously, corresponds to a single point on the- 
ses surfaces. Our calculation, shown in Fig. (5] provides the full amplitude -F'o(POiPjPs) 
in a two-dimensional domain. Computing this quantity is the main result of this work. 

4 Conclusion 

We have presented the solutions of the BS equation for the scattering states in Minkowski 
space. Results were limited to the spinless case, S-wave, and ladder kernel. The full 
off-mass-shell amplitude is found for the first time. Its on-mass shell values provide the 
phase shifts. They are in agreement with those obtained in Euclidean space by other 
methods. Relativistic phase shifts considerably differ, even at low energy, from the non- 
relativistic ones calculated by the Schrodinger equation with a Yukawa potential. Our 
results include the inelasticities which appear above the meson creation threshold. The 
knowledge of the off-mass-shell amplitude in Minkowski space is necessary to calculate 
the transition form factor and in the three-body BS-Faddeev equations. 
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